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Simulation  of  pores  sealing  during  homoepitaxy  on  Si(lll)  surface 

I.  G.  Neizvestny,  N.  L.  Shwartz,  Z Sh.  Yanovitskaja  and  A.  V.  Zverev 
Institute  of  Semiconductor  Physics  RAS  SB, 

1 3 pr.  Lavrenteva,  Novosibirsk,  630090,  Russia 

Abstract.  Simulation  of  homoepitaxy  process  on  porous  Si(  111)  using  3D  Monte  Carlo  model 
was  carried  out  to  give  estimation  of  necessary  dose  for  complete  pores  overgrowth  of  different 
sizes,  porosity  and  various  deposition  flux  intensities  for  given  temperature. 


Porous  silicon  (PS)  is  one  of  the  most  promising  materials  for  Si-based  emitting  devices 
that  can  be  useful  in  integrated  optoelectronics  [ , ].  Recently  an  approach  to  fabricate 
a compliant  substrate  for  growing  high-quality  heteroepitaxial  layers  was  proposed  [ ]. 
A double  layer  of  PS  obtained  from  heavily  doped  p-type  Si  was  used  as  the  substrate. 
A low-porosity  top  layer  and  a bottom  layer  with  high  porosity  act  as  a mechanical  damper 
between  heterostructure  and  a massive  support  and  allow  to  grow  GeSi  film  without  misfit 
and  threading  dislocation.  The  same  idea  was  used  when  growing  GaAs  layer  on  thin  Si 
layer  over  porous  Si  substrate  [ ].  Usually  there  are  only  circumstantial  evidences  that 
porous  substrate  is  overgrown  [ , ].  Number  of  deposited  layers  necessary  for  complete 
pores  healing  is  of  great  importance  for  highquality  heteroepytaxial  growth  on  porous 
substrates. 

3D  model  of  epitaxy  on  the  surfaces  of  diamond  like  crystals  [ ] was  applied  for  analysis 
of  epitaxy  on  (111)  porous  silicon  surfaces.  Our  lattice  model  can  take  into  account  atoms 
from  surface  region  up  to  hundred  atomic  layers  in  the  depth.  One  monolayer  could  contains 
up  to  105  atoms.  The  main  event  in  the  model  is  diffusion  hop,  with  atom  hopping  to  the 
empty  lattice  site  in  the  first,  second  or  third  coordination  spheres.  Adsorption  process 
is  considered  to  be  independent  on  diffusion.  Voids  and  vacancies  in  the  bulk,  multilevel 
steps  and  deep  pores  on  the  surface  are  possible.  The  abandonment  of  solid-on-solid  (SOS) 
principle  pennits  atoms  movement  along  the  walls  of  steps  and  pores. 

Simulation  of  homoepitaxy  process  on  porous  Si(  1 1 1)  was  carried  out  giving  estimation 
of  necessary  dose  for  complete  pores  overgrowth  of  different  sizes,  porosity  and  deposition 
rate.  Simulations  were  carried  out  in  the  following  ranges  of  substrate  parameters:  porosity 
13  < Ps  < 64%,  pores  diameter:  3.5-14  nm  and  rather  wide  ranges  of  deposition  rate, 
according  to  experimental  works  | , , ].  In  the  model  deposition  rate  is  introduced  through 
the  universal  parameter  n&a  ~ D/V,  where  D is  diffusion  coefficient  and  V-deposition 
rate.  Simulation  results  of  epitaxial  growth  on  porous  surface  with  100  x 100  atomic 
sites  (a.s.)  at  T — 1073  K [5,  ] for  different  surface  porosity,  pores  sizes  and  deposition 
rate  V = 0.004—40  nm/min)  are  presented  in  Fig.  1-4.  For  given  temperature  one  can 
estimate  deposition  rate  according  to  the  formulae  V — v ■ exp {—E&/kT)/n<M,  where  E d 
is  activation  energy  for  surface  diffusion  (eV),  V is  deposition  rate  (bilayer  per  second, 
BL/s),  v is  vibration  frequency  (1/s),  k is  Boltzmann’s  constant,  T is  substrate  temperature 
(K).  Taking  v — 1013  s"1,  £d  = 1.75  eV  [ ],  T = 1073  Kforndif  = 105  1/BL  one  obtains 
V = 0.6  BL/s. 
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Fig.  1.  Homoepitaxy  simulation  on  porous  Si(  111)  surface  at  T = 1073  K for  w^if  = 10  (high 
deposition  rate).  Pores  diameter  is  10  atomic  sites,  height  is  20  BL,  surface  porosity  Ps  = 26  %. 
(a)  initial  surface  top  view;  (b)  initial  cross-section  AA  along  [01 1 ] direction;  (c)  cross-section  AA 
after  1 .5  BL  deposition;  (d)  cross-section  AA  after  3. 1 BL  deposition. 


Fig.  2.  Dependence  of  minimum  number  of  deposited  bilayers  necessary  for  complete  pores 
overgrowth  Nc  on  porosity  Ps  for  small  pores  for  high  and  low  deposition  rates:  pores  diameter 
lOa.s.  T = 1073  K,  (1)  — ndif  = 10  (V  = 6000  BL/s),  (2)  — «dif  = 105  (V  = 0.6  BL/s). 


Plan  view  and  cross-section  along  [Oil]  direction  for  the  surface  with  small  pores  and 
surface  porosity  Ps  = 26%  before  and  after  deposition  with  high  deposition  rate  could  be 
seen  in  Fig.  I.  For  these  conditions  it  was  necessary  3.1  BL  of  silicon  to  be  pores  completely 
overgrown.  Figure  2 demonstrates  dependence  of  minimum  number  of  deposited  bilayers 
necessary  for  complete  pores  healing  ( Nc ) on  surface  porosity  Ps  for  low  and  very  high 
deposition  rates  and  small  pores.  For  high  deposition  rate  Nc  is  independent  of  porosity  and 
for  low  V Nc  increases  with  Ps.  This  dependence  is  clear;  for  high  deposition  rate  atoms 
reaching  pores  are  gathered  only  from  small  region  around  each  pore  limited  by  migration 
length.  For  porosity  less  than  65%  and  deposition  rates  under  examination  migration  length 
is  less  than  average  distance  between  pores.  For  low  deposition  rate  all  adatoms  on  the 
surface  are  involved  in  healing  process,  so  the  higher  porosity  the  less  numbers  of  atoms  are 
left  on  the  surface  (portion  of  atoms  falls  to  the  bottom  of  the  pores).  Figure  3 demonstrates 
dependence  of  Nc  on  ntn / for  three  different  pores  sizes.  For  large  pores  strong  dependence 
of  flux  intensity  is  observed.  Large  number  of  deposited  layers  Nc  is  necessary  for  complete 
pores  sealing  for  high  deposition  rate  (Nc  > 40  BL)  while  for  low  flux  Nc  is  rather  small 
and  has  only  weak  dependence  on  pores  diameter.  In  the  latter  case  Nc  is  determined 
by  deposited  dose  rather  than  diffusion,  so  dose  limited  regime  takes  place.  Figure  4 
illustrates  Nc  dependence  on  porosity  for  different  pores  sizes  (Fig.  4(a))  and  on  pores  size 
for  different  deposition  rates  (Fig.  4(b)).  Nc  independence  of  porosity  for  high  deposition 
rate  is  evident  from  curves  Fig.  4(a)  and  was  discussed  above.  Linear  dependence  of  Nc  on 
pore  diameter  for  low  deposition  rates  is  clear  from  curves  2,  3 while  curve  1 («dif  — 10, 
corresponds  to  high  deposition  rate)  demonstrates  superlinear  dependence.  Cross-sections 
along  [01 1 1 direction  before  (a)  and  after  deposition  (b-c)  on  porous  Si(l  1 1)  surface  with 


NT.13p 


135 


Fig.  3.  Dependence  of  Nc  on  «d if  for  different  sizes  of  pores.  T = 1073  K,  porosity  Ps  = 26%, 
pores  diameter:  (1 ) — 40  a.s.,  (2)  — 20  a.s.,  (3)  — 10  a.s. 


Fig.  4.  T = 1073  K (a)  Dependence  of  Nc  on  porosity  Ps  for  high  deposition  rate:  (naif  = 10)  for 
different  sizes  of  pores:  pores  diameter:  (1)  — 40  a.s.,  (2)  — 20  a.s.,  (3) — 10  a.s.;  (b)  Dependence 
of  Nc  on  pores  size  for  porosity  Ps  = 26%  and  different  growth  rate:  (1)  — Wdif  = 10,  (2)  — 
«dif  = 103,  (3)  — «dif  = 105. 


(a)  (b)  (c) 


Fig.  5.  Homoepitaxy  simulation  on  porous  Si(  111)  surface  at  T = 1073  K for  high  and  low 
deposition  rates,  surface  porosity  Ps  = 50%,  pores  diameter  40  a.s.,  height  20  BL.  (a)  initial  cross- 
section  along  [Oil]  direction;  (b)  cross-section  after  37.7  BL  deposition  for  high  deposition  rate 
(«dif  = 10);  (c)  cross-section  after  9.6  BL  deposition  for  low  deposition  rate  «dif  = 105. 


surface  porosi ty  Ps  = 50%  and  large  pores  are  presented  in  Fig.  5.  The  differences  in  sealing 
process  are  evident  from  these  figures.  For  high  deposition  rate  (Fig.  5(b))  overgrown  pores 
have  cone  shape  with  convex  up  bottom  and  Nc  is  found  to  be  rather  large.  In  this  case 
migration  length  is  small  and  adatoms  create  islands  as  between  pores  on  the  surface  as  on 
pores  bottom.  Atom  flux  in  our  simulation  is  perpendicular  to  the  surface.  Overhanging 
walls  arising  during  growth  process  are  beneficial  for  whiskers  formation  along  pores  walls 
near  the  bottom.  Quite  different  situation  one  has  for  low  deposition  rate  (Fig.  5(c)).  At 
the  initial  stage  of  growth  pores  shape  keeps  nearly  invariable,  only  diameter  of  each  pore 
is  decreasing  due  to  larger  migration  length  of  adatoms  mooving  along  pores  walls.  Final 
sizes  of  buried  pores  in  the  latter  case  are  noticeable  smaller. 
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